KSI)-TI)I{-65-274 


04 


lO 

vO 

t 

cc: 

Q 

H 

I 


U4 


a. 

8 

U4 

k-4 


l-H  I 

H  ‘‘ 
CO  • 
UJ  * 


CoD  RECORD  COPY 

RETURN  TO 

SCIENTIFIC  &  TECHNICAL  INFORMATION  DIVISION 
(ESTI),  BUILDING  1211 


esd  accession  list 

ESTI  Call  No._  /y  [  L/  ^  ^ 


Copy  No. 


Technical  Note 


1965-28 


B.  F.  LaPage 


Some  Properties 
of 

a  Scanned  Circular  Array 


7  July  1965 


Prepared  under  Electronic  Systems  Division  Contract  AF  19(628)-5167  by 


Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
Lexington,  Massachusetts 


h%1r,Liir 


The  work  reported  in  this  document  was  performed  at  Lincoln 
Laboratory,  a  center  for  research  operated  by  Massachusetts 
Institute  of  Technology,  with  the  support  of  the  U.S,  Air  Force 
under  Contract  AF  19(628)'5167. 


69 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 


SOME  PROPERTIES  OF  A  SCANNED  CIRCULAR  ARRAY 


TECHNICAL  NOTE  1965-28 

7  JULY  1965 


LEXINGTON  MASSACHUSETTS 


ABSTRACT 

The  impedance,  pattern  and  gain  properties  of  a  circular 
array  of  dipoles  (or  isotropic  sources)  is  discussed.  In  par¬ 
ticular  the  variation  of  these  properties  as  the  antenna  beam 
is  scanned  is  investigated.  Computed  results  for  typical  ar¬ 
rays  are  shown  and  the  meaning  of  the  results  are  discussed, 
and  further  areas  of  investigation  are  indicated. 


Accepted  for  the  Air  Force 
Stanley  J.  Wisniewski 
Lt  Colonel,  USAF 
Chief,  Lincoln  Laboratory  Office 


111 


f 


Some  Properties  of  a  Scanned  Circular  Array 


It  has  been  proposed  that  a  useful  satellite  antenna  design  could  be 
achieved  by  using  a  circular  array  of  monopoles  located  about  the  spin  axis  of 
a  short  cylinder  as  shown  in  Fig.  1.  The  spin  axis  of  the  satellite  is  oriented 
normal  to  the  orbital  plane  and  the  antenna  beam  is  **de-spun"  by  proper 
phasing  of  the  currents  on  the  array  elements.  The  array  elements  are  fed 
in  pairs  with  equal  amplitude  currents  from  separate  sources,  and  phased  to 
produce  the  desired  beam  positions.  A  retro-directive  scheme  could  also  be 
used,  but  in  any  event  some  basic  data  was  needed  on  the  behavior  of  scanned 
circular  arrays  and  this  was  the  motivation  of  this  study. 

There  is  a  fair  amount  of  literature  on  the  subject  of  circular  arrays'; 

1  2 

notably  Knudsen*s  exhaustive  study  and  more  recently  the  work  of  Neff, 

3  4 

Simpson  and  Mack.  Unfortimately,  these  studies  are  concerned  with  equi- 
phased  arrays,  parasitic  array  systems,  or  other  types  which  are  not  directly 
applicable  to  our  problem.  The  problem  of  the  input  impedance  variations  of 
an  element  during  scan  and  the  array  gain  variations  are  also  not  covered, 
and  these  topics  are  our  main  interest. 

5 

The  gross  parameters  of  the  array  were  chosen  from  what  was  con¬ 
sidered  to  be  a  reasonable  size  satellite  based  on  total  weight,  solar  cell 
output,  booster  capabilities,  etc;  and  are  typical,  at  least,  for  the  near  future. 

The  geometry  of  the  array  that  is  under  discussion  is  shown  in  Fig.  2. 
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The  N-elements  are  arranged  on  a  circle  of  radius  a,  located  in  the  x-y  plane 
and  phased  to  point  the  resultant  beam  at  an  angle  ~  90°).  Each  element 

is  positioned  by  the  angle  coordinate  (j)^  where 


(n-l) 


The  current  on  the  i  element  is  given  by 


j^i 

Ii=  e  ,  IIJ=  1  X 


where 


lb.  =~kacos(ct>  —  6  )  i 
o  n 


th 

Now  the  field  due  to  the  i  element  is 


E.  =  e 
1 


lb.  +  ka  cos((t>  —  d)  ) 

1  n  J 


or 


E.  =  e 
1 


j  ka  cos((t)  —  6  )  —  cos((b  —  cb  ) 

|_  n  o  n  J 


(1) 


(2) 


(3) 


The  total  field  is  obtained  by  summing  all  E.,  or 


N  jka 


E(<|>)  =  ^  e 
n=  1 


cos(<t)  — d)  )  —  cos(<|>  —  4)  ) 
n  o  n  J 


This  can  be  written  in  the  equivalent  form  as 


(4) 


2 


E(4,)=  y  e 
n='^l 


N  . 

j  u  sin  or 


n 


(5) 


<!>+  6 


O'  =  cf)  — 


n  n  2 


(j)  -  <|) 


u  =  2  ka  sin 


We  see  from  Eq.  (4)  that  when  (b  =  cb  ,  or  the  beam  points  to  an  element 

o  n 

position  angle,  E((j>)  does  not  depend  on  cj)  .  This  means  that  any  variatioils 


that  occur  in  E((|))  have  a  periodicity  of  the  inter-element  angular  spacing. 


and 


we  need  to  study  only  half  that  interval  to  obtain  all  the  scan  pattern  variations, 

i 

Now  Eq,  (5)  may  be  transformed  by  using  Bessel's  expansion: 


jZ  sin0 


UJ 

= 1 


(Z) 


n=  -CO 


(6) 


putting  the  result  in  a  form,  that  is  for  large  N,  easier  to  analyze.  That  is. 


E(<b)  J  (u)  +  2J  (u)  cos  —  (<t>  +  <|>  ),  N  even 
o  N  2  o 


E{<t))  =!  J  (u)  —  2jJ  (u)  sin  ^  (<t)  +  «|)  ),  N  odd 
o  N  2  o 


17) 

(8) 


u  =  2  ka  sin 


Now  we  show  that  for  N  S  6,  E(<|))  is  independent  of  N  to  a  large  degree.  At 


3 


the  null  of  E(()>)  we  have  J  (x)  =  0  for  x  =  2.  404  and  J-(2.  404)  =  .016, 

o  D 

JAZ.  404)  =  .  0033,  etc. 

6 

So  it  is  evident  that  terms  containing  J/{u)  are  negligible,  at  least  at  the 

D 

null  of  E {(})).  We  can  say  that  for  N  >  6  all  the  array  patterns  null  at  nearly 
the  same  angle,  given  by 

(b 

u  =  2kasin*^=  2.404 
o  2 

The  3  db  beamwidth  can  be  computed  similarly,  i.  e.  , 

(BW3)°=  4.in-‘  (i^)  . 


Also,  the  first  side  lobe  will  have  a  level  of  —7.  9  db  located  at  an  angle 


^  .  -1  ,3.  83, 

(si'e*  -  'tII’ 


We  can  now  give  the  array  plane  pattern  for  all  ring  arrays  with  N  ^  6,  as 
shown  in  Fig.  3.  For  N  =  4,  5  the  beamwidth  and  null  position  will  be  nearly 
the  same  as  in  Fig.  3  but  the  side  lobe  level  will  be  completely  different.  It 
should  be  noted  that  the  condition  u  ^  2ka  is  necessary  for  cj)  to  be  real,  and 
hence  the  pattern  to  be  in  the  visible  region.  The  pattern  of  the  array  in  the 
plane  normal  to  the  array  plane,  for  isotropic  sources,  can  be  analyzed  in  a 
similar  manner  to  obtain  the  result 


4 


N  even 


E(0,ct)  ) 

o 


J  (u)  +  J„(u) 
o  N 


TT  ,  ,  ,  X 


f  2  COS  N  +  <})  ) 
2  o 


(9) 


2j  sin  N  (—  +  <j)  ) 

L  O  J 


N  odd 


u  =  ka  (sin0  —  1 ) 


This  is  the  same  form  previously  given  except  for  a  change  in  the  argument 
of  J^(u).  One  might  ask  —  how  can  we  get  directivity  in  a  plane  orthogonal  to 
the  array  plane  for  isotropic  sources?  The  answer  is  seen  when  one  '^collapses 
the  point  sources  along  a  line  and  it  is  observed  that  we  have  an  end  fire  array 
with  tapered  amplitude  distribution.  This  is  known  to  produce  two-plane 
directivity;  the  well  known  Yagi  antenna  being  an  example  of  this  principle. 

Of  course,  if  array  elements  such  as  dipoles  are  considered,  then  the  array 
pattern  must  be  modified  to  take  into  account  the  element  pattern.  Figures 
4  and  5  show  the  array  3  db  beamwidth  vs.  array  diameter  (ka)  for  the 
*'in-plane”  and  '*out-of-plane**  cases.  We  can  now  estimate  the  gain  of  th^ 
antenna  array  by  using  the  well  known  expression 


G  ^ 


27,000 
^  o  ^  o 


where  0  ,  0  ,  are  the  principal  plane  half-power  beamwidths,  in  degrees. 

ill  id 

Using  the  data  in  Figs.  4  and  5  we  get  the  result  shown  in  Fig.  6  as  a  gain 
estimate.  Notice  that  the  gain  given  is  referenced  to  the  gain  of  an  element. 
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These  curves  tell  us  that  for  a  satellite  of  about  IX  diameter,  we  can  expect 
beamwidths  not  less  than  60®  wide  and  total  gains  of  10  -  12  db  over  isotropic. 
(This  assumes  that  monopoles  are  used  and  the  top  and  bottom  arrays  are 
spaced  0.  7X.  )  More  exact  gain  calculations  presented  later  show  that  the  gain 
estimates  of  Fig.  6  are  quite  good. 

Equation  (4)  has  been  used  to  machine  compute  a  series  of  patterns  as  a 
function  of  ka,  N,  cj),  aind  scan  angle  (j>^.  The  following  combinations  of  values 
were  used 

Ka  =  2,  4,  6,  8,  10  N=  4,  6,  7,  8,  10 

k  =  0,  .  1,  .  2,  .  3,  .  4,  .  5,  0  <  cl>  <  180%  Acj)  =  10® 

k  is  a  parameter  that  tells  us  where  the  beam  is  pointed,  i.  e.  , 

4)  =  k(f>  = 

A  range  of  0  <  k  <  .  5  scans  the  beam  one -half  an  inter -element  angular 
spacing  and  this  includes  all  the  desired  pattern  data  as  mentioned  before. 
Both  in-plane  and  out-of-plane  patterns  were  computed,  but  our  main  interest 
is  in  the  array  plane  patterns.  Some  of  this  data  is  given  in  Fig.  7  (a  through 
d).  The  patterns  for  N  >  6  show  the  general  characteristics  for  the  main 
beam  given  by  the  approximate  formulas,  but  as  expected  the  side  and  back 
lobe  levels  are  considerably  different.  The  patterns  for  N  =  4,  Ka  =  2.  28 
(Fig.  8)  for  example,  show  a  back  lobe  for  some  scan  angles  whose  intensity 
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equals  that  of  the  main  beam.  With  a  beamwidth  of  58 compared  to  56"^  from 
Fig.  3,  we  see  that  the  beamwidth  approximation  is  still  good  for  N  =  4.  For 
the  detailed  study  of  an  array  pattern,  the  more  exact  formula  (Eq.  5)  is 
necessary,  but  the  approximations  are  still  generally  quite  useful. 

We  now  turn  our  attention  to  the  problem  of  the  variation  of  the  input  im¬ 
pedance  of  each  dipole  element  during  scan  and  also  a  calculation  of  the  gain 
of  the  array.  Figure  9  shows  the  dipole  array  geometry.  Since  the  impedance 
of  a  monopole  over  a  infinite  groimd  plane  is  one -half  that  of  a  dipole  the  re¬ 
sults  of  the  calculations  here  are  easily  modified  to  include  the  monopoly  case. 
The  input  impedance  of  an  array  element  is  given  by  the  circuit  equations 


V  I  I 

z  = -^  =  z  +  z  z 

1  11  12  In 


2^9  =  2;  -i+  z  + - +  z 

2  21  22  2n 

2  2 


Z  = 
n 


Z 

nn 


(10) 


Let  the  power  input  to  element  s  be  and  the  far-field  intensity  generated  be 

E  .  We  assume  E  =  kl  where  I  is  the  element  current, 
s  s  s  s 
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Then 


Z  I'sl' W 

s 

and  for  a  single  element 

P  =  |I  R  (Z  ) 
o  I  o  I  e  oo 


The  field  intensities  are  related  by 


E  = 
o 


k|I 


o 


Since  the  fields  are  "phased*'  to  produce  a  maximum,  the  power  gain  is  defined 


as 


G  = 


—  for  P  =  P_ 
2  o  T 


and  we  have 


Now  in  our  case  I  =  e  ,  since  III  =  1. 

S  I  S  I 
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Then 


G  = 


R  (Z  ) 
e  oo 


y  R  (z ) 

l_i  e  s 


Since 


R  (Z  )  =  y  -flZ  .  I  cos  (or  .  +  (j>.  —  <|)  n- 
e  s  sjl  sj  j  s  J 

j 


where 


a  . 
SJ 


tan 


X  . 
SJ 

R 

SJ 


Z 

SJ 


R  .  +  jX  . 

SJ  SJ 
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phase  of  the  current  in  j 


th 


element. 


The  final  form  of  the  gain  expression  is: 


G  (over  an  element)  = 


N^R  (Z  ) 
_ e  oo 

N 


We  also  have  the  important  relations: 


Z  =  Z  (identical  elements)  , 

oo  ss 

Z  =  Z 


mn 


nm 


(reciprocity  of  mutual  impedance)  , 


(t)j  =  ka  cos  —  (p  —  i+  l)(0<p<  —  scan  factor) 

Now  Eq,  (11)  reduces  to  the  form  G  =  N  if  all  Z  =0  for  m  ^  n,  as  it  should. 

mn 

The  formula  also  agrees  exactly  with  some  results  given  in  Kraus^  for  two 
elements  spaced  one-half  wavelength,  with  in  phase  and  out-of-phase  excita¬ 
tions.  In  order  to  provide  more  flexibility,  a  separate  computer  subroutine 

is  used  to  compute  the  mutual  impedances,  given  the  length  of  an  element  L/\ 

7 

and  the  spacing  d/X.  The  induced  EMF  method  is  used,  which  assumes  sinus¬ 
oidal  antenna  currents  and  infinitely  thin  dipoles.  The  mutual  impedance  of 

g 

parallel  spaced  dipoles  can  be  expressed  by  using  the  Ci  and  Si  functions, 
but  for  machine  computation  more  effort  is  expended  in  calculating  these 
functions  than  in  a  direct  calculation  of  the  mutual  impedance  itself.  The  input 
parameters  to  the  main  computer  program  are  as  follows: 

Parameter  Values  Used 

N  -  number  of  elements  4,  5,  6 

ka  -  ring  size  2.  28,  2.  50,  2.  75,  3.  0 

L/\  -  dipole  length,  wavelengths  .  475,  .  500,  .  525 

The  program  computes  the  following: 

1.  All  mutual  impedances  Z  .. 

s  J 

2.  All  self  impedances  Z 

oo 

3.  Input  impedance  of  each  element  vs.  scan  angle. 
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4.  Gain  of  the  array  vs.  scan  angle. 

A  typical  calculation  result  is  shown  in  Fig.  10.  The  gain  variation  is 
seen  to  be  very  small.  However,  the  input  impedance  variation  is  quite  large. 
In  fact,  some  combinations  of  parameters  even  lead  to  negative  real  values  of 
input  impedance.  This  indicates  that  more  power  is  flowing  towards  the  gen¬ 
erator  for  that  element,  than  away  from  it.  An  operating  condition  that  one 
would  certainly  wish  to  avoid.  To  show  the  variations  in  a  more  readily  im- 
derstandable  format,  we  have  normalized  the  dipole  input  impedance  values 
to  100  ohms  and  plotted  them  on  a  Smith  Chart.  This  is  equivalent  to  normal¬ 
izing  a  monopole  input  impedance  to  a  50  ohm  transmission  line,  a  practical 
consideration  in  any  actual  system.  Figures  1 1  (a  through  m)  shows  some 
typical  data  along  with  the  gain  variations  and  self  impedance  values.  The 
impedance  variation  causes  the  VSWR  to  vary  between  2  and  3.  If  we  are 
willing  to  give  up  the  continuous  beam  scanning  in  favor  of  a  step- scanned 
system,  we  can  avoid  the  ends  of  the  impedance  variations  curve  and  reduce 
the  variation  to  about  a  VSWR  of  1.7.  This  step- scanning  procedure  causes 
discrete  overlapping  beams  and  a  subsequent  loss  in  gain  under  some  condi¬ 
tions.  For  the  N  =  6  cajse,  for  instance,  the  stepping  would  be  in  60®  incre¬ 
ments  and  would  produce  a  3.  0  db  ripple  in  the  scanned  pattern.  Another 
possible  way  to  minimize  the  impedance  variation  is  to  use  isolators  in  each 
element  feed  line.  A  VSWR  of  2.  0  is  effectively  reduced  to  1. 23  for  10  db 
isolation  loss  and  for  20  db  loss  the  effective  VSWR  is  1.07.  Recent  advances 
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in  the  design  of  coaxial  isolators  has  made  available  \mits  that  have  20  db  re¬ 
verse  loss  and  0.  7  db  forward  loss  and  weigh  less  than  one  po\md.  These 
values  are  typical  of  commercial  units  and  have  not  been  optimized,  especially 
as  regards  the  weight  of  the  unit,  for  satellite  use. 

It  is  the  author's  opinion  that  an  analysis  of  the  properties  of  a  ring  array 
of  dipoles  much  more  rigorous  than  that  presented  here  is  of  little  immediate 
value.  The  reason  for  this  is  that  the  mathematical  model  of  monopoles  on  an 
infinite  ground  plane  bears  a  questionable  resemblance  to  monopoles  arrayed 
on  the  top  and  bottom  of  a  small  cylinder.  Until  an  exact  analysis  of  the  actu¬ 
al  geometry  is  made  and  then  checked  with  measured  data,  the  correlation 
between  model  and  real  world  is  not  really  known.  It  is  felt  that  a  survey 
type  look  at  a  large  number  of  possible  satellite  antenna  schemes  is  more  de¬ 
sirable,  so  that  the  most  promising  of  these  may  receive  more  intensive 
scrutiny  from  the  analytical  and  the  experimental  viewpoint. 

Some  computations  concerned  with  the  allowable  tolerance  on  the  excita¬ 
tion  currents  necessary  to  keep  gain  variations  within  some  arbitrary  bound 
have  been  made  by  the  author,  and  the  results  are  given  in  Lincoln  Laboratory 
Technical  Memorandum  61L-0016. 

Mr.  Robert  Peck  (Group  61)  programmed  the  IBM- 7094  for  the  computa¬ 
tions  given  in  this  report. 
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Fig.  Z.  Array  geometry. 
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Zira/\  =  ka 


I.OX  I.5X  2.0X  2.5X 

DIAMETER 

Fig.  4.  Ring  array  beamwidth  vs.  diameter. 
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HALF -POWER  BEAMWIDTH  (deg) 


I.OX  I.5X  2.0X  2.5X 

DIAMETER 

Fig.  5.  Ring  array  beamwidth  vs.  diameter. 
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Fig.  6.  Ring  array  gain  vs.  diameter. 
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RELATIVE  POWER  (db) 


Fig.  7-a.  Pattern  of  circular  array 
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RELATIVE  POWER  (db) 


Fig.  7-c.  Pattern  of  circular  array. 


RELATIVE  POWER  (db) 


Fig.  7-d.  Pattern  of  circular  array. 
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Fig.  9. 


Dipole  array  geometry. 
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Fig.  10.  Input  impedcince  vs.  scan  angle. 
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IMPEDANCE  COORDINATES 
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Fig.  11-b. 


Ring  array  impedance. 
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2  =  63.1  +  j  10.4 


G  =  6.335  ±0.015  db 


Fig.  11-c.  Ring  array  impedance. 
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IMPEDANCE  COORDINATES 


.  11-d. 


Fig 


Ring  array  impedance 


IMPEDANCE  COORDINATES 


73.1  + i42.5 


13-61-3647 

N  =  6 
ka  =  2.50 
L  =  0.500 


G=  6.715  ±  0.095  db 


Fig.  11-e. 


Ring  array  impedaince. 
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IMPEDANCE  COORDINATES 


84.7  +j  75.5 


13-61-36481 


Fig.  11-f. 


Ring  array  impedance. 
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IMPEDANCE  COORDINATES 
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IMPEDANCE  COORDINATES 


Fig.  11-h. 


Ring  array  impedance. 
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IMPEDANCE  COORDINATES 


13-61-3651 


G  =  6.05±0db 


N=5 

ka  =  2.28 
L  =  0.525 


Fig.  11-i. 


Ring  array  impedance. 
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IMPEDANCE  COORDINATES 


Fig.  11-j. 


Ring  array  impedance. 
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IMPEDANCE  COORDINATES 


Fig.  11-k. 


Ring  array  impedance. 
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63.1+ j  10.4 


G  =  5.195  ±0.345db 


Fig.  11-1.  Ring  array  impedajice. 
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Z„  =  73.l  +  j42.5 


3-61-3655 


N  =  4 
ka=2.28 
L  =  0.500 


G  =  5.19  ±  0.34  db 


Fig.  11-m. 


Ring  array  impedance. 
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